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Abstract 

J/^-nuclear bound state energies are calculated for a range of nuclei by solving the Proca 
(Klein-Gordon) equation. Critical input for the calculations, namely the medium-modified D 
and D* meson masses, as well as the nucleon density distributions in nuclei, are obtained from 
the quark-meson coupling model. The attractive potential originates from the D and D* meson 
loops in the J/^f self-energy in nuclear medium. It appears that J/^-nuclear bound states should 
produce a clear experimental signature provided that the J/^> meson is produced in recoilless 
kinematics. 
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I. INTRODUCTION 



With the 12 GeV upgrade of the CEBAF accelerator at the Jefferson Lab in the USA and 
with the construction of the FAIR facility in Germany, we expect tremendous progress in 
understanding the properties of charmonia and charmed mesons in nuclear medium. These 
new facilities will be able to produce low-momenta charmonia and charmed mesons such 
as J/\l/, ip(2S), D and D* in an atomic nucleus. Because the targets are nuclei and the 
nuclear fermi-momentum is available, it may also be possible to produce these mesons at 
subthreshold energies. While at JLab charmed hadrons will be produced by scattering 
electrons off nuclei, at FAIR they will be produced by the annihilation of antiprotons on 
nuclei. One of the major challenges is to find appropriate kinematical conditions to produce 
these mesons essentially at rest, or with small momentum relative to the nucleus. One of 
the most exciting experimental efforts may be to search for the J/^-nuclear bound states 
amongst many other possible interesting experiments in these facilities. The discovery of 
such bound states would provide evidence for a negative mass shift of the J/\l/ meson, and 
a possible role of the QCD color van der Waals forces [1] in nuclei. Ref. [2] presents a recent 
review of the properties of charmonium states and compiles a complete list of references for 
theoretical studies concerning a great variety of physics issues related to these states. 

There is a relatively long history for the studies of charmonium binding in nuclei. The 
original suggestion [1] that multiple gluon QCD van der Waals forces would be capable of 
binding a charmonium state led to an estimate of a binding energy as large as 400 MeV in 
an A = 9 nucleus. However, using the same approach but taking into account the nucleon 
density distributions in the nucleus, Ref. [3] found a maximum of 30 MeV binding energy 
in a large nucleus. Along the same line, including the nuclear density distributions in nuclei 
r) c - nuclear bound state energies were estimated in Ref. [4]. Based on Ref. [5], which showed 
that the mass shift of charmonium in nuclear matter can be expressed in terms of the usual 
second-order Stark effect arising from the chromo-electric polarizability of the nucleon, the 
authors of Ref. [6] obtained a 10 MeV binding for J/\l> in nuclear matter, in the limit of 
infinitely heavy charm quark mass. On the other hand, for the excited charmonium states, a 
much larger binding energy was obtained, e.g. 700 MeV for the ip'{2S) state, an admittedly 
untrustworthy number. Following the same procedure, but keeping the charm quark mass 
finite and using realistic charmonium bound-state wave-functions, Ref. [7] found 8 MeV 
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binding energy for J/\I/ in nuclear matter, but still over 100 MeV binding for the charmonium 
excited states. While a larger value for the QCD Stark effect is expected for excited states 
(because of the larger sizes of these states), the large values for the binding energies found in 
the literature for these states may be an overestimate. A possible source of this overestimate 
is the breakdown of the multipole expansion for the larger-sized charmonium states. 

There are some other studies of charmonium interactions with nuclear matter, in par- 
ticular involving the J/^f meson. QCD sum rules studies estimated a J/$ mass decrease 
in nuclear matter ranging from 4 MeV to 7 MeV [8-10], while an estimate based on color 
polarizability [11] gave larger than 21 MeV. Since the J/\I> and nucleons have no quarks in 
common, the quark interchange or the effective meson exchange potential should be negli- 
gible to first order in elastic scattering [1], and multi-gluon exchange should be dominant. 
Furthermore, there is no Pauli blocking even at the quark level. Thus, if the J/\l/-nuclear 
bound states are formed, the signal for these states will be sharp and show a clear, narrow 
peak in the energy dependence of the cross section, as will be explained later. This situation 
has a tremendous advantage compared to the cases of the lighter vector mesons [12]. 

Previously, we have studied [13] the J/\I/ mass shift (scalar potential) in nuclear matter 
based on an effective Lagrangian approach, including the D and D* meson loops in the J/ty 
self-energy. This is a color-singlet mechanism at the hadronic level, and may be compared 
to the color-octet mechanism of multi-gluon exchange, or QCD van der Waals forces. In this 
work, we extend our previous study made in nuclear matter [13] to finite nuclei, and com- 
pute the J/\l/-nuclear bound state energies by explicitly solving the Proca (Klein-Gordon) 
equation for the J/\I/ meson produced nearly at rest. Furthermore, the structure of heavy 
nuclei, as well as the medium modification of the D and D* masses, are explicitly included 
based on the quark- meson coupling (QMC) model [14]. 

A first estimate for the mass shifts (scalar potentials) for the J/^f meson (and \&(3686) 
and \I/(3770)) in nuclear medium including the effects of D meson loop in the J/ty self- 
energy - see Fig. 1 - was made in Ref . [7] . Employing a gauged effective Lagrangian for the 
coupling of D mesons to the charmonia, the mass shifts were found to be positive for J /iff 
and ■0(3770), and negative for -0(3660) at normal nuclear matter density p . These results 
were obtained for density-dependent D and D masses that decrease linearly with density, 
such that at p they are shifted by 50 MeV. The loop integral in the self-energy (Fig. 1) 
is divergent and was regularized using form-factors derived from the 3 P decay model with 
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FIG. 1: DD-loop contribution to the J/\E' self-energy. We include also DD* and D*D* contribu- 
tions. 

quark- model wave functions for ip and D. The positive mass shift is at first sight puzzling, 
since even with a 50 MeV reduction of the D masses, the intermediate state is still above 
threshold for the decay of J/ty into a DD pair and so a second-order contribution should 
be negative. However, as we have shown in Ref. [13], this is a result of the interplay of 
the form factor used and the gauged nature of the interaction used in Ref. [7]. We have 
estimated [13] the mass shift (scalar potential) of the J/ty meson including the DD, DD*, 
D*D and D*D* meson loops in the self-energy using non-gauged effective Lagrangians. The 
density dependence of the D and D* masses is included by an explicit calculation [15] using 
the quark- meson coupling (QMC) model [16]. The QMC model is a quark-based model 
for nuclear structure which has been very successful in explaining the origin of many-body 
or density dependent effective forces [17] and hence describing nuclear matter saturation 
properties. It has also been used to predict a great variety of changes of hadron properties 
in nuclear medium, including the properties of hypernuclei [18]. A review of the basic 
ingredients of the model and a summary of results and predictions can be found in Ref. [19]. 

II. J/\& POTENTIAL IN INFINITE NUCLEAR MATTER 

We briefly review how the J/^f scalar potential is calculated in nuclear matter [13]. We use 
the effective Lagrangian densities at the hadronic level for the vertices, J/^/DD, J/^/DD* 
and J /^D*D* (in the following we denote by ip the field representing J/^f): 
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C^dd = ig^DDV [D (d^D) - (d^D) D] , (1) 
C^dd- = e a ^ u (d a ^) \(d»D*») D + D (d^D* v )} , (2) 

£<pd*d* = ig^D* D * {V [(d„D w ) d; - D* u {d^DD] 

+ [(d,r) d: - r fao:)] d^ + d^ ir - (<vn m}- (3) 

The Lagrangian densities above are obtained on the basis of SU(4) invariance for the cou- 
plings among the pseudo-scalar and vector mesons [20, 21]. We use the values for the 
coupling constants, g^DD = g^D*D* = g^DD* = 7.64, where the former two values were 
obtained using the vector meson dominance model [20]. The form of the J/^D*D* coupling 
respects Yang's theorem [22, 23] and so for spin one J/\I> the virtual D* and D* must not 
be simultaneously transverse. The difference with the gauged Lagrangian of Ref. [7] for the 
i\)DD vertex amounts to adding to Eq. (1) a term 2g"i DD ijj ll 'ijj IM DD. It may be of interest to 
investigate in future whether the meson loops investigated here might play some role in the 
understanding of the well known J/ip — >■ prr puzzle [24]. 

Only the scalar part contributes to the J/\& self-energy, and the scalar potential for the 
J/ty meson is the difference of the in- medium, m^, and vacuum, m^, masses of J/^f, 

Am = m* 4> - rrty, (4) 

with the masses obtained from 

m% = K) 2 + £(A; 2 = ml) . (5) 

Here is the bare mass and S(/c 2 ) is the total J/ty self-energy obtained from the sum 
of the contributions from the DD, DD* and D*D* loops. The in-medium mass, m^, is 
obtained likewise, with the self-energy calculated with medium-modified D and D* meson 
masses. 

In the calculation, we use phenomenological form factors to regularize the self-energy 
loop integrals following a similar calculation for the p self-energy [25]: 

For the vertices ^DD, ^DD* and ^D*D*, the form factors used are respectively, F DD (q 2 ) = 
u 2 D (q 2 ), F DD *(q 2 ) = u D (q 2 ) u D *(q 2 ), and F D * D *(q 2 ) = u 2 D ,(q 2 ) with A D and A D * being the 



cutoff masses, and the common values, Ad = Ad* are used. For the calculation in finite 
nuclei, we chose two values within the range expected to be reasonable, Ad = Ad* = 
1500 and 2000 MeV to estimate the model dependence. The scalar potential, (m^ — m^), 
calculated via Eq. (5) in nuclear matter [13], is shown in Fig. 2 as a function of nuclear 
matter density. 
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FIG. 2: The scalar potential arising from DD, DD* and D*D* loops as a function of nuclear 
matter density for different values of the cutoff mass Ad = A^* [13]. 



III. J/* BOUND STATES IN FINITE NUCLEI 

In this section we discuss the case that J/i/j is placed in a finite nucleus. The nucleon den- 
sity distributions for 12 C, 16 0, 40 Ca, 90 Zr and 208 Pb are obtained using the QMC results for 
finite nuclei [14]. For 4 He, we use the parametrization for the density distribution obtained 
in Ref. [26]. Then, using a local density approximation we calculate the J/\l/ potentials in 
nuclei. Note that also the in-medium masses of D and D* mesons, which are necessary to 
estimate the J/\l/ self-energy consistently, are calculated in the QMC model. We emphasize 
that the same quark-meson coupling constants are used between the applied meson mean 
fields and the light quarks in the nucleon and those in the D and D* mesons [15], where 
these coupling constants are calibrated by the nuclear matter saturation properties [19]. 

As examples, we show the J/\l/ potentials calculated for 4 He and 208 Pb nuclei in Fig. 3 
for the two values of the cutoff mass, 1500 MeV and 2000 MeV. The maximum depth of the 
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J/ty potential ranges roughly from 18 MeV to 24 MeV. 
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FIG. 3: J potential in a 4 He nucleus (left panel) and in a 208 Pb nucleus (right panel) for two 
values of the cutoff mass in the form factors, 1500 and 2000 MeV. 



Using the J potentials in nuclei obtained in this manner, we next calculate the J/^f- 
nuclear bound state energies. We follow the procedure applied in the earlier work [27] for 
the rj- and cu-nuclear bound states. In this study, we focus entirely on the situation where 
the J/ty meson is produced nearly at rest (recoilless kinematics in experiments). Then, it 
should be a very good approximation to neglect the possible energy difference between the 
longitudinal and transverse components [28] of the J/\l/ wave function (4>%,) as was assumed 
for the u meson case [28]. After imposing the Lorentz condition, <9^0^, = 0, to solve the 
Proca equation, aside from a possible width, becomes equivalent to solving the Klein-Gordon 
equation with the reduced mass of the system /i in vacuum [29, 30]: 



[V 2 + EI-/J 2 - 2^{m%{r) - m 9 )] fa(r) = 0, 



(7) 



where is the total energy of the J/ty meson, and /i = m^M/(m^ + M) with m$ (M) 
being the vacuum mass of the J/\I/ meson (nucleus). Since in free space the width of J/ty 
meson is ~ 93 keV [31], we can ignore this tiny natural width in the following. When the 
J/ty meson is produced nearly at rest, its dissociation process via J/^f + N — > + D is 
forbidden in the nucleus (the threshold energy in free space is about 115 MeV above). This 
is because the same number of light quarks three participates in the initial and final states 
and hence, the effects of the partial restoration of chiral symmetry which reduces mostly 
the amount of the light quark condensates, would affect a similar total mass reduction for 



7 



the initial and final states [15, 32]. Thus, the relative energy (~ 115 MeV) to the threshold 
would not be modified significantly. 

We also note that once the J/^f meson is bound in the nucleus, the total energy of 
the system is below threshold for nucleon knock-out and the whole system is stable. The 
exception to this is the process J/ty + N — > r] c + N which is exothermic. Ref. [33] has 
provided an estimate of this cross section from which we deduce a width for the bound J/^f 
of order 0.8 MeV (nuclear matter at ps = Po/2). We therefore expect that the bound J/^, 
while not being completely stable under the strong interaction, should be narrow enough 
to be clearly observed. It would be worth while to investigate this further. Then, provided 
that experiments can be performed to produce the J/\l/ meson in recoilless kinematics, the 
J/\I/ meson is expected to be captured by the nucleus into one of the bound states, which 
has no strong interaction originated width. This is a completely different and advantageous 
situation compared to that of the rj and uj meson cases [27]. Thus, in the end, we may simply 
solve the Klein-Gordon equation with the reduced mass Eq. (7), without worrying about the 
width, under the situation we consider now. The calculation is carried out in momentum 
space by the method developed in Ref. [29, 30]. Bound state energies for the J/^f meson, 
obtained solving the Klein-Gordon equation Eq. (7), are listed in Table I. 

The results in Table I show that the J/\l/ are expected to form nuclear bound states in 
all the nuclei considered. This is insensitive to the values of cutoff mass used in the form 
factors. It is very interesting to note that one can expect to find a J/ , I'- 4 He bound state. 
It will be possible to search for this state as well as the states in a 208 Pb nucleus at JLab 
after the 12 GeV upgrade. In addition, one can expect quite rich spectra for medium and 
heavy mass nuclei. Of course, the main issue is to produce the J/*f> meson with nearly 
stopped kinematics, or nearly zero momentum relative to the nucleus. Since the present 
results imply that many nuclei should form J/\I/-nuclear bound states, it may be possible to 
find such kinematics by careful selection of the beam and target nuclei. 

IV. SUMMARY AND DISCUSSION 

In this study, we have calculated the J/\l/-nuclear potentials using a local density approx- 
imation, with the inclusion of the DD, DD* and D*D* meson loops in the J/\I/ self-energy. 
The nuclear density distributions, as well as the in-medium D and D* meson masses, are 
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TABLE I: Bound state energies calculated for the J/^> meson, solving the Klein-Gordon equation 
with the reduced mass Eq. (7). Under the situation that the J/^f meson is produced in recoilless 
kinematics, the width due to the strong interactions is all set to zero, as well as its tiny natural 
width of ~ 93 keV [31] in free space. (See also the text.) 





A DtD . = 1500 MeV 


A D D , = 2000 MeV 




E (MeV) 


E (MeV) 


iHe 


Is 


-4.19 


-5.74 




Is 

ID 


-9.33 


-11.21 




1 n 


-2.58 


-3.94 


16 Q 


1 S 

ID 


-11.23 


-13.26 




1 n 


-5.11 


-6.81 




Id 


-14.96 


-17.24 




IP 

XT 


-10.81 


-12.92 




Id 


-6.29 


-8.21 




2s 


-5.63 


-7.48 


|°Zr 


Is 


-16.38 


-18.69 




lp 


-13.84 


-16.07 




Id 


-10.92 


-13.06 




2s 


-10.11 


-12.22 


| 08 Pb 


Is 


-16.83 


-19.10 




lp 


-15.36 


-17.59 




Id 


-13.61 


-15.81 




2s 


-13.07 


-15.26 



consistently obtained employing the quark-meson coupling model. In the model, all the cou- 
pling constants between the applied meson mean fields and the light quarks in the nucleon, 
as well as those in the D and D* mesons, are all equal and calibrated by the nuclear matter 
saturation properties. Using the J/\l/-nuclear potentials obtained in this manner, we have 
solved the Klein-Gordon equation which is reduced from the Proca equation, and obtained 
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the J/\I>-nuclear bound state energies. For this, we have been able to set the strong inter- 
action width of the J /iff meson to be zero (and neglect its tiny natural width of ~ 93 keV 
in free space). Combined with the generally advocated color-octet gluon-based attraction, 
or QCD color van der Waals forces, we expect that the J /iff meson will form nuclear bound 
states and that the signal for the formation should be experimentally very clear, provided 
that the J /iff meson is produced in recoilless kinematics. 
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